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Abstract 



Results are reported from a search for the effects of contact interactions using events 
with a high-mass, oppositely-charged muon pair. The events are collected in proton- 
proton collisions at a/s = 7 TeV using the Compact Muon Solenoid detector at the 
Large Hadron Collider. The data sample corresponds to an integrated luminosity of 
5.3 fb^^. The observed dimuon mass spectrum is consistent with that expected from 
the standard model. The data are interpreted in the context of a quark- and muon- 
compositeness model with a left-handed isoscalar current and an energy scale param- 
eter A. The 95% confidence level lower limit on A is 9.5 TeV under the assumption of 
destructive interference between the standard model and contact-interaction ampli- 
tudes. For constructive interference, the limit is 13.1 TeV. These limits are comparable 
to the most stringent ones reported to date. 
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1 Introduction 

The existence of three famiUes of quarks and leptons might be explained if these particles are 
composed of more fundamental constituents. In order to confine the constituents (often re- 
ferred to as "preons" [|1J |2i) and to account for the properties of quarks and leptons, a new 
strong gauge interaction, metacolor, is introduced. Below a given interaction energy scale A, 
the effect of the metacolor interaction is to bind the preons into metacolor-singlet states. For 
parton-parton center-of-mass energy less than A, the metacolor force will manifest itself in the 
form of a flavor-diagonal contact interaction (CI) [S'/l]. In the case where both quarks and lep- 
tons share common constituents, the Lagrangian density for a CI leading to dimuon final states 
can be written as 



C 



(1) 



qi ={go/ ^^){VLL{qL7^'qL){}lL7liFL) 

+ ?/rl(ur7^uk)(/Jl7hj"l) 

+ ?/RL(dR7^'dR)(/iL7FFL) 
+ ?/Rr(ur7''ur) iflRjpijiR) 

+ '/RR(dR7^dR)(f/R7j,;^R)}, 

where qi = (u, d)L is a left-handed quark doublet, ur and dR are right-handed quark singlets, 
and }iL and }Ir are left- and right-handed muons. By convention, gQ/4in = 1. The parameter 
A characterizes the compositeness energy scale. The parameters rjij allow for differences in 
magnitude and phase among the individual terms. Lower limits on A are set separately for 
each term with rjij taken, by convention, to have a magnitude of one. 

The dimuons from the subprocesses for standard model (SM) Drell-Yan (DY) (Si production 
and from CI production can have the same helicity state. In this case, the scattering amplitudes 
are summed, resulting in an interference term in the cross section for pp — )■ X + as 
illustrated schematically in Fig.[l] 

The differential cross section corresponding to the combination of a single term in Eq. [l] with 
DY production can be written as 



I 



c 



dM 



dM 



mi 



(2) 



where M^i is the invariant dimuon mass, X is due to interference, and C is purely due to the 
CI. Note that rjij = +1 corresponds to destructive interference and rjij = — 1 to constructive 
interference. The processes contributing to the cross section in Eq. (|2| are denoted collectively 
by "CI/DY". The difference da^^'^"^ / dM^, - da^^/dM^ is the signal we are searching for in 
this paper. 





Figure 1: Schematic representation of the addition of DY (left) and CI (right) amplitudes, for 
common helicity states, contributing to the total cross section for pp — > X + ji^ji^ . 



2 



2 Predictions of the left-left isoscalar model 



The contact interaction model used for this analysis is the left-left isoscalar model (LLIM) IH, 
which corresponds to a left-handed current interaction described by the first term of £q/ in 
Eq. (1). The LLIM is the conventional benchmark for CI in the dilepton channel. For this 
analysis, all initial-state quarks are assumed to be composite. 

Previous searches for CI in the dijet and dilepton channels have all resulted in limits on the com- 
positeness scale A. Searches have been reported from experiments at LEP ||6l- [T0l , HERA IITTl 
[m, the Tevatron lUMSl, and recently from the ATLAS iSH22I and CMS Il23ti25l experiments 
at the LHC. The best limits in the LLIM dimuon channel are A > 9.6 TeV for destructive inter- 
ference and A > 12.9 TeV for constructive interference, at the 95% confidence level (CL) [22J. 

In this paper, we report a search for CI in the dilepton channel produced in pp collisions at 
i/s = 7 TeV using the Compact Muon Solenoid (CMS) detector at the Large Hadron Collider 
(LHC). The data sample corresponds to an integrated luminosity of 5.3 fb^^. 

2 Predictions of the left-left isoscalar model 

The basic features of the LLIM dimuon mass spectra are demonstrated with a generator-level 
simulation using PYTHIA ||26| , with appropriate kinematic selection criteria that approximate 
the acceptance of the detector. Figures |2]ja) and |2]jb) show the LLIM dimuon mass spectra for 
different values of A for destructive and constructive interference, respectively. The curves 
illustrate that with increasing mass the CI leads to a less steeply falling yield relative to DY 
production, with the effect steadily increasing with decreasing A. For a given value of A, the 
event yield is seen to be larger for constructive interference compared to the destructive case, 
with the relative difference increasing with A. 

For the results presented in this paper, the analysis is limited to a dimuon mass range from 




M (GeV/c^) (GeV/c^) 



Figure 2: Simulated dimuon mass spectra using the left-left isoscalar model for different values 
of A for (a) destructive interference and (b) constructive interference. The events are generated 
using the PYTHIA Monte Carlo program with kinematic selection requirements that approxi- 
mate the acceptance of the detector. As A increases, the effects of the CI are reduced, and the 
event yield approaches that for DY production. The model predictions are shown over the full 
mass range, although the model is not valid as M;,^,c^ approaches A. The integrated luminosity 
corresponds to 63 fb^^. 



3 



200 to 2000GeV/c2. The lower mass is sufficiently above the Z-peak so that a deviation from 
DY production would be observable. The highest dimuon mass observed is between 1300 and 
1400 GeV/ and, for the values of A where the limits are set, less than one event is expected 
for dimuon masses above 2000GeV/c^. In order to limit the mass range in which the detector 
acceptance has to be evaluated, we therefore choose an upper mass cutoff of 2000GeV/c^. To 
optimize the limit on A, a minimum mass M^l^ is varied between the lower and upper mass 
values, as described in Section |6] 

3 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam- 
eter, providing a magnetic field of 3.8 T. Within the field volume are a silicon pixel and strip 
tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass/ scintillator hadron 
calorimeter. Muons are measured in gas-ionization detectors embedded in the steel flux-return 
yoke. Extensive forward calorimetry complements the coverage provided by the barrel and 
endcap detectors. A detailed description of the CMS detector can be found in Ref. [27J. 

The tracker and muon detector are important sub-systems for this measurement. The tracker 
measures charged particle trajectories within the range |//| < 2.5, where pseudorapidity tj = 
— ln[tan(0/2)], and polar angle d is measured from the beam axis. The tracker provides a 
transverse momentum (px) resolution of about 1% at a few tens ofGeV/c to 10% at several 
hundred GeV/c [28 1, where px is the component of momentum in the plane perpendicular to 
the beam axis. Tracker elements include about 1400 silicon pixel modules, located close to 
the beamline, and about 15 000 silicon microstrip modules, which surround the pixel system. 
Tracker detectors are arranged in both barrel and endcap geometries. The muon detector com- 
prises a combination of drift tubes and resistive plate chambers in the barrel region and a com- 
bination of cathode strip chambers and resistive plate chambers in the endcap regions. Muons 
can be reconstructed in the range < 2.4. 

For the trigger path used in this analysis, the first level (LI) selects events with a muon candi- 
date based on a sub-set of information from the muon detector. The trigger muon is required 
to have \r]\ < 2.1 and px above a threshold that was raised to 40GeV/c by the end of the data- 
taking period. This cut has little effect on the acceptance for muon pairs with masses above 
200GeV/c^. The small effect is included in the simulation. The high level trigger (HLT) refines 
the LI selection using the full information from both the tracker and muon systems. 

4 Event selection criteria 

This analysis uses the same event selection as the search for new heavy resonances in the 
dimuon channel, discussed in Ref. Il29ll . Each muon track is required to have a signal ("hit") 
in at least one pixel layer, hits in at least nine strip layers, and hits in at least two muon de- 
tector stations. Both muons are required to have px > 45GeV/c. To reduce the cosmic ray 
background, the transverse impact parameter of the muon with respect to the beamspot is re- 
quired to be less than 0.2 cm. In order to suppress muons coming from hadronic decays, a 
tracker-based isolation requirement is imposed such that the sum of px of all tracks, excluding 
the muon and within a cone surrounding the muon, is less than 10% of the px of the muon. The 
cone is defined by the condition AK = y^C^'/I^'+T^^)^ — ^-3, where (p is the azimuthal angle 
of a track, and the differences A;/ and A^ are determined with respect to the muon's direction. 

The two muons are required to have opposite charge and to be consistent with originating 
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5 Simulation of SIV! and CI dimuon production 



Table 1: Description of event samples with detector simulation. The cross section a and inte- 
grated luminosity L are given for each sample generated. 





Process 


Generator 


# of Events 


t7(pb) 


L(pb-^) 


Order 


7/'v* - 


-4 1I1J M,,„ > ^?ncpv/r^ 


PYTHIA 




7.90x10" 


6.91x10'' 


LO 


Z/y* - 


-> }i}i,M,,,, > 200GeV/c2 


PYTHIA 


5.50x10* 


9.70x10^1 


5.67x10* 


LO 


Z/7* - 


-> }i}i,Mfi^ > 500GeV/c2 


PYTHIA 


5.50x10" 


2.70x10^2 


2.04x101^ 


LO 


Z/7* - 


}i}i,Mmi > SOOGeV/c^ 


PYTHIA 


5.50x10* 


3.10x10^3 


1.77x10^ 


LO 


Z/7* - 


> }i}i, Mfif, > lOOOGeV/c^ 


PYTHIA 


5.50x10" 


9.70x10-4 


5.67x10^ 


LO 




Z/7* TT 


PYTHIA 


2.03x10^ 


1.30 xlO^ 


1.56x10^ 


LO 




tt 


MadGraph 


2.40x10^ 


1.57 xlO^ 


1.54x10^ 


NLO 




tw 


POWHEG 


7.95 xl05 


7.90x10° 


1.01 xlO^ 


NLO 




tw 


POWHEG 


8.02 xlO^ 


7.90x10° 


1.02x10^ 


NLO 




WW 


PYTHIA 


4.23x10'' 


4.30 xlO^ 


9.83x10* 


LO 




WZ 


PYTHIA 


4.27x10'' 


1.80 xlO^ 


2.37x10^ 


LO 




zz 


PYTHIA 


4.19x10'' 


5.90x10° 


7.10x10^ 


LO 




W + jets 


MadGraph 


2.43x10^ 


3.10x10" 


7.82 xlO^ 


NLO 


multi-jet, }i (pt > 15 GeV/c) 


PYTHIA 


1.08x10^ 


8.47x10* 


1.28x10^ 


LO 



from a common vertex. To suppress cosmic ray muons that are in time with the collision event, 
the angle between the two muons must be smaller than n — 0.02 radians. At least one of the 
reconstructed muons must be matched (within AR < 0.2 and Apj/ pj < 1) to the HLT muon 
candidate. 

If an event has more than two muons passing the above requirements, then the two highest-px 
muons are selected, and the event is retained only if these muons are oppositely charged. Only 
three such events are observed with selected dimuon mass above 200GeV/c^, and in all three 
cases, the dimuon mass is less than 300 GeV/c^. Thus, events with multiple dimuon candidates 
play essentially no role in the analysis. 



5 Simulation of SM and CI dimuon production 

This section describes the method used to simulate the mass distribution from the CI/DY pro- 
cess of Eq. (|2|, including the leading-order (LO) contributions from DY and CI amplitudes, their 
interference, the effects of next-to-leading-order (NLO) QCD and QED corrections, and the re- 
sponse of the detector. The predicted number of CI/DY events is the product of the generated 
number of CI/DY events, a QCD K-factor, a QED K-factor, and a factor denoted as "acceptance 
times migration" (A x M). The factor A x M is determined from the detector simulation of DY 



events, as explained below in Section 5.2 The simulation of background due to non-DY SM 
processes is also described. 

5.1 Event samples with detector simulation 

A summary of the event samples used for simulation of the detector response to various 
physics processes is presented in Table [T] The event generators used are PYTHIA, with the 
CTEQ6.6M implementation [30J of parton distribution functions (PDF), POWHEG ll3Tti33l . and 
MadGraph5 IMI- The detector simulation is based on Geant4 [35J. 



5.2 Detector acceptance times mass migration 

To simplify the analysis, we use the detector simulation for DY events to determine the detector 
response for CI/DY events, which have a behavior similar to that for DY events for the large 
values of A of interest in this analysis. For a given value of MJJJ]", the product of acceptance 



5.2 Detector acceptance times mass migration 
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Figure 3: Acceptance times migration, A x M, versus MJJJ]". Corresponding values and uncer- 
tainties are given in Table |2] The error bars indicate statistical uncertainties based on simulation 
of the DY process. The systematic uncertainty is 3% as explained in the text. The increase of 
A X M at lower mass is due to the increase in acceptance, while at higher mass, it is dominated 
by the growth in mass resolution. Since the cross section falls steeply with mass, events tend to 
migrate from lower to higher mass over a range determined by the mass resolution. 

times migration (A x M) is given by the ratio of the number of DY events reconstructed with 
mass above M'^ to the number of DY events generated with mass above MJJ^'". Some of the 
reconstructed events have been generated with mass below M|Jj™ because of the smearing due 
to the mass reconstruction, which has a resolution of 6.5% at masses around 1000 GeV/c^, rising 
to 12% at 2000GeV/c2. The dependence of A x M on M^^j" is plotted in Fig. [s] and values are 
given in Table|2] The increase of A x M at lower mass is due to the increase in acceptance, while 
at higher mass, it is dominated by the growth in mass resolution. Since the cross section falls 
steeply with mass, events tend to migrate from lower to higher mass over a range determined 
by the mass resolution. 

To validate that the A x M factor based on DY production is applicable to CI/DY production, 
we compare event yields predicted using the A x M factor with those predicted using a simu- 
lation of CI/DY production. The study is performed for the cases of constructive interference 
with A = 5 and lOTeV, which represent a wide range of possible CI/DY cross sections. The 
results differ by at most 3%, consistent with the statistical precision of the study. The systematic 
uncertainty in A x M is conservatively assigned this value. 

5.2.1 Event pileup 

During the course of the 2011 data taking period, the luminosity increased with time, resulting 
in an increasing "event pileup", the occurrence of multiple pp interactions recorded by the de- 
tector as a single event. The dependence of reconstruction efficiency on event pileup is studied 
by weighting simulated events so that the distribution of the number of reconstructed primary 
vertices per event matches that in data. The reconstruction efficiency is found to be insensitive 
to the variations in event pileup encountered during the data taking period. 
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5 Simulation of SIV! and CI dimuon production 



Table 2: Multiplicative factors used in the prediction of the expected number of events from the 
CI/DY process. The uncertainties shown are statistical. The systematic uncertainty is 3% for 
A X M and 3% for the QCD K-factor, as explained in the text. The uncertainty in the QED K- 
factor is dominated by the systematic uncertainty that is assigned as the size of the correction, 
I (QED K-f actor) — 1 1, to allow for systematic uncertainty in the generator. 



M;j);"(GeV/c2) 


AxM 


QCD K-factor 


QED K-factor 


200 


0.80 ± 0.01 


1.303 ±0.005 


1.01 


300 


0.82 ± 0.01 


1.308 ±0.005 


0.99 


400 


0.83 ± 0.01 


1.299 ±0.005 


0.97 


500 


0.86 ±0.02 


1.305 ±0.005 


0.95 


600 


0.86 ± 0.01 


1.299 ±0.005 


0.94 


700 


0.87 ±0.01 


1.298 ±0.005 


0.92 


800 


0.88 ± 0.01 


1.288 ±0.005 


0.91 


900 


0.89 ± 0.01 


1.280 ±0.004 


0.90 


1000 


0.89 ± 0.01 


1.278 ±0.004 


0.89 


1100 


0.89 ± 0.01 


1.275 ±0.004 


0.88 


1200 


0.91 ±0.01 


1.268 ±0.004 


0.88 


1300 


0.92 ±0.01 


1.262 ±0.004 


0.87 


1400 


0.94 ±0.01 


1.260 ±0.004 


0.87 


1500 


0.97 ±0.01 


1.261 ±0.004 


0.86 



5.3 Higher-order strong and electromagnetic corrections 

Since we use the leading-order generator pythia to simulate the CI/DY production, we must 
determine a QCD K-factor which takes into account higher-order initial-state diagrams. Under 
the assumption that the QCD K-factor is the same for DY and CI/DY events, we determine the 
QCD K-factor as the ratio of DY events generated using the next-to-leading-order generator 
MC@NLO IMl to those generated using pythia. The MC@NLO generator is used with the 
same PDF set as used with PYTHIA. The resulting QCD K-factor as a function of M|JJf is given 
in Table |2] The large sizes of the simulated event samples result in statistical uncertainties less 
than 0.5%. The systematic uncertainty is assigned the value 3%, the size of the correction II37I 
between next-to-next-to-leading order (NNLO) and NLO DY cross sections. For SM processes 
other than DY production, the QCD K-factor is found, independent of dimuon mass, from 
the ratio of the cross section determined using MC@NLO to the cross section determined from 
PYTHIA. 

The effect of higher-order electromagnetic processes on CI/DY production is quantified by a 
mass-dependent QED K-factor determined using the HORACE generator Il38l . The values of 
the QED K-factor, as a function of M'^^, are given in Table |2] The systematic uncertainty is 
assigned as the size of the correction, | (QED K-factor) — 1 1, since the effect of higher-order QED 
corrections on the new physics of CI is unknown. 

5.4 Non-DY SM bacl<grounds 

Using the samples of simulated events listed in Table 111 event yields are predicted for various 
non-DY SM background processes, as shown in Table |3] The yields are given as a function of 
M|^J" and they are scaled to the integrated luminosity of the data, 5.28 ± 0.12 fb^^ |39|. For 
comparison, the expected yields are also shown for DY events. The relevant backgrounds, 
in decreasing order of importance, are tt, diboson (WW/WZ/ZZ), W (including W+jets and 
tW), and Z — > TT production. The background from multi-jet events is studied using both the 
simulation sample listed in Table [T] and control samples from data, as reported in Ref. P29II. 
The results of either method indicate that no multi-jet background events are expected for 
M™" > 200GeV/c2. For MJ]lf > lOOOGeV/c^ the fractional statistical uncertainty in the non- 



5.5 Predicted event yields 
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Table 3: Expected event yields for DY and non-DY SM backgrounds. The uncertainties shown 
are statistical. A systematic uncertainty of 2.2% arises from the determination of integrated 
luminosity |,39J. 



(GeV/c2) 


DY 


tt 


Diboson 


W+Jets & tW 


Z^- TT 


Sum Non-DY 


200 


3630 ± 18 


454 ±3 


123.0 ±2 


47.90 ±1.35 


6.96 ±4.14 


632.3 ±5.9 


300 


870.6 ± 8.8 


104 ±2 


38.6 ±1.2 


12.82 ±0.70 





155.9 ±2.1 


400 


301.6 ±5.1 


26.0 ±0.8 


12.7 ±0.7 


3.32 ±0.35 





42.0 ±1.1 


500 


123.8 ±3.3 


8.19 ±0.46 


5.07 ±0.41 


1.02 ±0.20 





14.3 ±0.6 


600 


55.31 ±0.19 


2.92 ± 0.27 


2.42 ± 0.28 


0.29 ±0.11 





5.63 ± 0.41 


700 


27.35 ± 0.13 


1.12 ±0.17 


0.86 ±0.16 


0.07 ±0.05 





2.06 ± 0.24 


800 


14.23 ±0.10 


0.34 ± 0.09 


0.51 ±0.12 


0.07 ±0.05 





0.92 ±0.16 


900 


7.72 ± 0.07 


0.05 ±0.03 


0.25 ± 0.08 


0.07 ±0.05 





0.36 ±0.10 


1000 


4.32 ± 0.05 


0.05 ± 0.03 


0.10 ±0.05 


0.07 ±0.05 





0.21 ±0.08 


1100 


2.46 ± 0.04 


0.05 ±0.03 


0.09 ± 0.05 


0.07 ±0.05 





0.20 ±0.08 


1200 


1.48 ± 0.03 





0.01 ±0.01 


0.07 ±0.05 





0.08 ±0.05 


1300 


0.91 ±0.02 





0.01 ±0.01 


0.07 ±0.05 





0.08 ±0.05 


1400 


0.56 ±0.02 





0.01 ±0.01 


0.07 ±0.05 





0.08 ±0.05 


1500 


0.33 ± 0.02 








0.07 ±0.05 





0.07 ±0.05 



DY background is large, but the absolute yield is much smaller than that for DY background. 

5.5 Predicted event yields 

Using the methods described above, the sum of the event yields for the CI/DY process and the 
non-DY SM backgrounds, for the integrated luminosity of the data sample, are predicted as a 
function of M^^^ and A. The predicted event yields for destructive and constructive interfer- 
ence are given in Tables |4] and |5] 

For destructive interference, there is a region of the M^^-A parameter space where the pre- 
dicted number of events is less than for SM production. This "reduced-yield" region is indi- 
cated in Table ji] The region of parameter space, MJJJf > 600GeV/c^ and A < 12TeV, where 
our expected limit is most stringent (see Fig. |5ja)), lies outside the reduced-yield region. For 
constructive interference, the predicted number of events is always larger than for SM produc- 
tion. 

6 Expected and observed lower limits on A 

6.1 Dimuon mass distribution from data 

The observed numbers of events versus MjJ^J^ are given in Table Wl The observed distribution 
of M^t is plotted in Fig. |4] along with the expected distributions irom the SM and for CI/DY 
plus non-DY SM processes, for three illustrative values of A. The data are consistent with the 
predictions from the SM, dominated by DY production. 

6.2 Limit-setting procedure 

Since the data are consistent with the SM, we set lower limits on A in the context of the LLIM. 
The expected and observed 95% CL lower limits on A are determined using the CLg modified- 
frequentist procedure described in Ii40.,41j , taking the profile likelihood ratio as a test statistic 
||42| . The expected mean number of events for a signal from CI is the difference of the number 
of CI/DY events expected for a given A, and the number of DY events. The expected mean 
number of background events is the sum of events from the DY process and the non-DY SM 
backgrounds. The observed and expected numbers of events are given in Tables |4] and |5] 
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6 Expected and observed lower limits on A 



Table 4: Observed and expected number of events for illustrative values of M^^. The expected 
yields are shown for SM production and for the sum of CI/DY production (for destructive 
interference and for a given A) and non-DY SM backgrounds. For each column of MJJ^'", the 
expected yield for CI/DY + non-DY SM production that is just above that expected for SM 
production is highlighted in red. Entries above the highlighted ones correspond to values of 
A for which the expected yield is less than that for SM production, because of the destructive 



interference term in the cross section. As discussed in Section 6.3 the best expected limit is 
obtained for MJJJJ" = llOOGeV/c^. For this choice, the expected event yield is highlighted in 
gray that corresponds to the value of A closest to the observed 95% CL lower limit on A of 
9.5 TeV (9.7TeV expected). 
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141 


57 


28 


14 
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8 


3 
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60.9 


29.4 


15.2 


8.1 


4.5 


2.7 


1.6 


1.0 


0.6 


0.4 


A (TeV) MC 
























18 


134.2 


58.0 


27.9 


14.3 


7.7 


4.3 


2.6 


1.5 


1.0 


0.6 


0.4 


17 


134.5 


57.9 


27.7 


14.4 


7.8 


4.4 


2.6 I 




1.0 


0.6 


0.4 


16 


134.9 


58.0 


27.8 


14.5 


7.8 


4.5 1 




1.6 


1.0 


0.7 


0.5 


15 


135.6 


58.3 


28.1 


14.7 


8.0 } 




2.9 


1.7 


1.1 


0.8 


0.5 


14 


133.7 


58.3 


28.3 


15.0 


8.4 


5.0 


3.1 


1.9 


1.3 


0.9 


0.6 


13 


134.1 


59.3 


29.1 


15.7 


8.9 


5.4 


3.5 


2.2 


1.5 


1.0 


0.7 


12 


138.6 


60.1 


30.2 


16.7 


9.8 
11.2 


6.1 


4.1 


2.7 


1.9 


1.3 


0.9 


11 


135.7 


62.5 


32.1 


18.4 


7.3 


5.0 


3.5 


2.4 


1.7 


1.2 


10 


141.1 


66.7 


35.7 


21.2 


13.6 


9.2 


6.6 


4.6 


3.3 


2.4 


1.7 


9 


148.5 


73.8 


42.4 


27.1 


18.3 


13.1 


9.5 


6.9 


5.0 


3.7 


2.6 


8 


164.7 


88.1 


54.4 


36.8 


26.2 


19.3 


14.3 


10.6 


7.8 


5.8 


4.1 


7 


198.1 


117.5 


79.4 


57.6 


43.3 


31.6 


24.0 


17.5 


13.1 


9.0 


6.2 


6 


278.1 


182.3 


131.7 


100.1 


7(,.7 


57.9 


45.1 


33.0 


22.6 


16.9 


11.3 


5 


469.2 


338.7 


261.7 


204.4 


158.6 


123.2 


96.7 


74.6 


56.8 


41.5 


29.2 


4 


1025 


784.1 


620.1 


494.2 


384.3 


302.6 


232.8 


174.6 


127.2 


94.5 


68.5 


3 


3199 


2517 


2012 


1599 


1242 


975.7 


744.7 


575.4 


437.7 


320.1 


231.4 



Systematic uncertainties in the predicted signal and background event yields are estimated 
from a variety of sources and included as nuisance parameters in the limit-setting procedure. 
Significant sources of systematic uncertainty are given in Table |6] The uncertainty in the in- 
tegrated luminosity is described in Ref. [39]. The uncertainty in the CI/DY acceptance is ex- 



plained in Section 5.2 The uncertainties in the prediction of backgrounds depend on the value 
of MJ]^. These uncertainties are given in Table ro^ for the values of MJJJ™ chosen for limits on 
A with destructive and constructive interference rThe PDF uncertainty in the expected yield of 
DY events is evaluated using the PDF4LHC procedure [43 J. The uncertainties in the QED and 



QCD K-factors are explained in Sec. 5.3 The uncertainty from non-DY backgrounds is due to 
the statistical uncertainty associated with the simulated event samples. The systematic uncer- 
tainties which decrease the limit on A by the largest amounts are the uncertainties on the PDF 
and QED K-factor. When both these uncertainties are set to zero, the limit for destructive in- 
terference is increased by 0.4% and the limit for constructive interference is increased by 3.0%. 
Thus, the systematic uncertainties degrade the limits by only small amounts. 

We considered possible systematic uncertainties in modeling the detector response by com- 
paring kinematic distributions between data and simulation of DY and non-DY SM processes. 
There are no differences in these distributions that could lead to significant systematic uncer- 
tainties through their effect on selection efficiency and mass resolution. 



6.3 Results for limits on A 
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Table 5: Observed and expected number of events as in Table |4| Here CI/DY predictions are 
for constructive interference. Shown with gray highlighting is the expected event yield corre- 
sponding to the value of A closest to the observed 95% CL lower limit on A of 13.1 TeV (12.9 TeV 
expected) for M^}^ selected to be 800GeV/c^. 
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3.1 


2.1 


1.4 


14 
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2.4 


11 
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5.5 
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41.6 
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11.4 


8.1 
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Table 6: Systematic uncertainties affecting the limit on A, evaluated for the values of MJJ|{" that 
provide the best expected limits for constructive and destructive interference. 





Uncertainty (%) 


Source 


const. 


destr. 


integrated luminosity 


2.2 


2.2 


acceptance times migration {A x M) 


3.0 


3.0 


PDF 


13.0 


16.0 


QED K-factor 


9.0 


11.8 


QCD K-factor 


3.0 


3.0 


DY MC statistics 


1.2 


1.6 


non-DY backgrounds 


1.1 


2.9 



6.3 Results for limits on A 

The observed and expected lower limits on A at 95% CL as a function of M^'" for destructive 
and constructive interference are shown in Figs.j5|a) and|5|b). The value of M^™, chosen to 
maximize the expected sensitivity, is llOOGeV/c^or destructive interference, and 800GeV/c^ 
for constructive interference. The observed (expected) limit is 9.5 TeV (9.7 TeV) for destructive 
interference and 13.1 TeV (12.9 TeV) for constructive interference. The variations in the observed 
limits lie almost entirely within the \-a (standard deviation) uncertainty bands in the expected 
limits, consistent with statistical fluctuations. The number of expected events corresponding to 
the observed limits on A are highlighted in Tables |4] and |5] 



7 Summary 

The CMS detector is used to measure the invariant mass distribution of pairs produced in 
pp collisions at a center-of-mass energy of 7 TeV, based on an integrated luminosity of 5.3 fb^^ . 
The invariant mass distribution in the range 200 to 2000GeV/c^ is found to be consistent with 
standard model sources of dimuons, which are dominated by Drell-Yan production. The data 
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7 Summary 



CMS, ^ = 7 TeV, 5.3 fb" 




M,, (GeV/c^) 



Figure 4: Observed spectrum of M^^ and predictions for SM and CI/DY plus non-DY SM 
production. Predictions are shown for three illustrative values of A, for constructive and de- 
structive interference. The error bars for data are 68% Poisson confidence intervals. 

are interpreted in the context of a quark- and muon-compositeness model with a left-handed 
isoscalar current and an energy scale parameter A. The 95% confidence level lower limit on A 
is 9.5 TeV under the assumption of destructive interference between the standard model and 
contact-interaction amplitudes. For constructive interference, the limit is 13.1 TeV. These limits 
are comparable to the most stringent ones reported to date. 
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Figure 5: Observed and expected limits as a function of MJJJ)" for (a) destructive interference 
and (b) constructive interference. The value of MJJJf^, chosen to maximize the expected sensitiv- 
ity, is 1100 GeV/c-^ for destructive interference, and 800 GeV/c^ for constructive interference. The 
observed (expected) limit is 9.5 TeV (9.7 TeV) for destructive interference and 13.1 TeV (12.9 TeV) 
for constructive interference. The observed limit at the value chosen for M™" is indicated with 
a red plus sign. The variations in the observed limits lie almost entirely within the l-cr bands, 
consistent with statistical fluctuations. 
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